Magnetoresistance due to edge spin accumulation 
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Because of spin-orbit interaction, an electrical current is accompanied by a spin current resulting 
in spin accumulation near the sample edges. Due again to spin-orbit interaction this causes a small 
decrease of the sample resistance. An applied magnetic field will destroy the edge spin polarization 
leading to a positive magnetoresistance. This effect provides means to study spin accumulation by 
electrical measurements. The origin and the general properties of the phenomenological equations 
describing coupling between charge and spin currents are also discussed. 



It was predicted a long time ago [J 0| that because 
of spin-orbit interaction electrical and spin currents are 
interconnected: an electrical current produces a trans- 
verse spin current and vice versa. In recent years this 
has become a subject of considerable interest. 

The purpose of this Letter is twofold. First, another 
way of understanding this interconnection will be pre- 
sented and some general properties of the resulting phe- 
nomenological equations will be discussed. Second, a new 
magnetoresistance effect will be considered, which allows 
to study the current-induced spin accumulation near the 
sample edges by purely electric measurements. 

The transport phenomena related to coupling of the 
spin and charge currents can be described phenomenolog- 
icaly in the following simple way. Let q be the electron 
flow density and let be its conventional expression 
not accounting for spin-orbit interaction: 

q (0) = -finE - DVn, (1) 

where /i and D are the usual electron mobility and dif- 
fusion coefficient, connected by the Einstein relation, E 
is the electric field, and n is the electron concentration. 
The electric current density is j = —eq, where e is the 
absolute value of the electron charge. 

Let qij be the spin polarization current density tensor 
(the flow of the j component of the spin polarization in 
the direction i). It should be understood that for polar- 
ized electrons the spin current may exist even in the ab- 
sence of spin-orbit interaction, simply because spins are 
carried by electron flow. We denote the corresponding 
quantity as q^K Then, similar to Eq. (1), we have 

fl—Wi-D^ (2) 

where P is the vector of electron spin polarization den- 
sity. If there are other sources for currents, like for ex- 
ample a temperature gradient, the corresponding terms 
should be included in Eqs. (1) and (2). 

We have departed from the conventional definitions [l|, 
Q by introducing the vector of spin polarization density P 
and the spin ■polarization current qij . This allows to avoid 
numerous factors 1/2 and 2 in the formulas to follow. 
One can return to the traditional notations by putting 
P = 2S, where S is the spin density, and by replacing 
Qij by Qy/2 to obtain the true spin current density. 



Spin-orbit interaction couples the two currents. For a 
material with inversion symmetry Q we have: 

(0) , (0) / Q x 

Qi = qi +7£ijkq) k ', (3) 

qij = qfk ~ 7%fe% , (4) 

where is the unit antisymmetric tensor and 7 is a di- 
mensionless coupling constant proportional to the spin- 
orbit interaction, it is assumed that 7 << 1. The differ- 
ence in signs in Eqs. (3) and (4) is consistent with the 
Onsager relations and is due to the different properties 
of q and qij with respect to time inversion [3] ■ 

Explicit phenomenological expressions for the two cur- 
rents follow from Eqs. (l)-(4): 

j/e = fmE + DX7n + (3E A P + 5 curl P, (5) 

dP dn 
= -fiEiPj - D—L + eijk{(3nE k + S-^-). (6) 

Here 

= 7/ i, S = jD, (7) 

so that the coefficients j3 and 5, similar to fj, and D, satisfy 
the Einstein relation. 

Eqs. (5) and (6) should be complemented by the equa- 
tion for the spin polarization vector: 

^ + H + (S!AF)j + & = 0, (8) 

where the vector Jl is directed along the applied magnetic 
field, SI being the spin precession frequency and r s is 
the spin relaxation time. In Eqs. (6), (7) we ignore 
the action of magnetic field on the particle dynamics. 
This is justified if lo c t << 1, where oj c is the cyclotron 
frequency and r is the momentum relaxation time. Since 
normally r s >> r, it is possible to have both f2r s >> 1 
and uj c t << 1 in a certain interval of magnetic fields. It 
is also assumed that the equilibrium spin polarization in 
the applied magnetic field is negligible. 

While Eqs. (5)-(8) are written for a three-dimensional 
sample, they are equally applicable to the 2D case, with 
obvious modifications: the electric field, space gradients, 
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and all currents (but not the spin polarization vector) 
should have components in the 2D plane only. 

In the equilibrium situation all currents should obvi- 
ously vanish. If an inhomogeneous magnetic field ex- 
ists, the equilibrium spin polarization will be space- 
dependent, however this by itself should produce neither 
spin, nor charge currents. To assure this, an additional 
counter-term should be introduced into the right-hand 
side of Eq. (2), proportional to dBj/dxi, which takes 
care of the force acting on the electron with a given spin 
in an inhomogeneous magnetic field B(r) (see [5]). Cor- 
responding terms will appear in Eqs. (5), (6). We ignore 
these terms assuming that B is homogeneous. 

Equations (5)-(8), which appeared for the first time in 
Refs. describe all the physical consequences of spin- 

charge current coupling [6 L The term (3E A P describes 
the anomalous Hall effect 7] , where the spin polarization 
plays the role of the magnetic field. 

The term 5 curl P describes an electrical current in- 
duced by an inhomogeneous spin density (now referred 
to as the Inverse Spin Hall Effect). A way to measure 
this current under the conditions of optical spin orienta- 
tion was proposed in The circularly polarized excit- 
ing light is absorbed in a thin layer near the surface of 
the sample. As a consequence, the photo-created electron 
spin density is inhomogeneous, however curlP = 0, since 
both P and its gradient are perpendicular to the surface. 
By applying a magnetic field parallel to the surface one 
can create a parallel component of P, thus inducing a 
non-zero curl P and the corresponding surface electric 
current (or voltage). This effect was found experimen- 
tally for the first time by Bakun et al @. 

The term (3neijkEk (and its diffusive counterpart 
Seijkdn/dxk) in Eq. (6), describes what is now called the 
Spin Hall Effect: an electrical current induces a trans- 
verse spin current, resulting in spin accumulation near 
the sample boundaries [l], [2] . This phenomenon was ob- 
served experimentally only in recent years [Tol . [Tl| and 
has attracted widespread interest. 

It should be stressed that all these phenomena are 
closely related and have their common origin in the cou- 
pling between spin and charge currents given by Eqs. (3) 
and (4). Any mechanism that produces the anomalous 
Hall effect will also lead to the spin Hall effect and vice 
versa. It is remarkable that there is a single dimension- 
less parameter, 7, that governs the resulting physics. The 
calculation of this parameter should be the objective of 
a microscopic theory. For the case, when the coupling is 
due to spin asymmetry in electron scattering, this was 
done in Ref. [2J, where (3 and 5 were expressed through 
the scattering amplitude. In this case 7 depends only on 
the form of the scattering potential, the electron energy, 
and the strength of spin-orbit interaction. 

An "intrinsic" mechanism of the spin Hall effect, re- 
lated only to spin band splitting, was proposed for bulk 
holes in the valence band [T2|. The value of 7 is on the 
order of (fcp^) -1 (kp is the Fermi wavevector, I is the 
mean free path) and generally depends on the details of 




FIG. 1: The geometry of the proposed experiment with a 2D 
sample. The direction of the magnetic field is of no impor- 
tance, so long as it lies in the 2D plane. For simplicity it is 
assumed to be parallel to the current. 



the scattering mechanism. The current consensus is that 
the intrinsic mechanism may exist for any type of spin 
band splitting, except if it is linear in k [13[ . 

Note that the J = 3/2 holes may not be described 
by the simple Eqs. (5)-(8), because for higher spins the 
number of coupled macroscopic quantities increases com- 
pared to spin 1/2 particles. The mutual transformation 
of spin and charge currents for holes, due to scattering, 
was studied in [14j . Also, even in the absence of spin- 
orbit interaction, holes are still particles with internal 
angular momentum L = 1 and the splitting into light 
and heavy holes still exists. Thus, for the case of holes 
the spin-orbit interaction is not of primary importance. 

We now discuss a new related phenomenon: a mag- 
netoresistance due specifically to spin accumulation near 
the sample edges. Since the accumulation occurs on the 
scale of the spin diffusion length L s — \JDt s (the "spin 
layer" [l|), the proposed effect depends on the sample 
size, L, and becomes negligible when L » L s . 

Whithin the spin layer the z component of spin po- 
larization changes in the direction perpendicular to the 
sample boundary (the y direction). Thus curl_P ^ 0, and 
according to Eq. (5) a correction to the electric current 
should exist. As we will see, this correction is positive, 
i.e. it leads to a slight decrease of the sample resistance 
compared to the (hypothetical) case when spin-orbit in- 
teraction is absent. By applying a magnetic field in the 
xy plane, we can destroy the spin polarization (the Hanle 
effect) and thus observe a positive magnetoresistance on 
a field scale corresponding to Slr s ~ 1. One might say 
that this is a manifestation of combined direct and in- 
verse spin Hall effects, and the Hanle effect. 

We will consider a 2D sample (see Fig. 1), a similar 
effect will exist also for a thin wire. The advantage of the 
2D case is that the small effect considered here will not 
be masked by the normal magnetoresistance, because the 
magnetic field parallel to the 2D plane acts on the spins 
only, but not on the electron orbital motion. Since the 
spin polarization is proportional to the electric field, we 
discard nonlinear in E terms proportional to EP. 



P(B) 




10 12 14 16 18 20 



n-c. 



FIG. 2: Normalized magnetoresistance as a function of the 
parameter flr s for different sample widths. 1 - Eq. (18) for 
A » 1, 2 - A = 1.5, 3 - A = 0.8, 4 - A = 0.5, 5 - A = 0.3 



For the geometry of Fig. 1, from Eq. (5) we obtain: 



HP 

j = e(imE + 5^), 
dy 



so that the total current is 

r-L/2 



j{y)dy = J + AI, 

-L/2 



(9) 



(10) 



where 

I = ennEL, AI = e5(P z (L/2) - P z {—L/2)). (11) 

The correction to the current, AI, is proportional to the 
difference in spin polarization at the opposite edges of 
the sample. Eq. (6) yealds: 



dy 



D ^ < 12 » 



In the steady state Eq. (8) gives: 



d 2 P P 

dy 2 y t s 



d 2 P P 

D±^ + np z = ^. 

dy 2 t s 



(13) 



These equations should be solved with the boundary con- 
ditions at y = ±L/2: 



dP z _ (3nE 
~dy~~~' 



dy 



(14) 



corresponding to vanishing spin currents q yz and q yy at 
the sample edges. 

A straightforward calculation gives the result: 



where i?o is the uncorrected sample resistance, AR is the 
field-dependent negative correction due to spin accumu- 
lation, 



k = (1 — ix) 



1/2 



x = Qt s 



A = L/(2L S ), 



and L s = {Dt s ) x / 2 is the spin diffusion length. 

Thus AR is proportional to the square of the dimen- 
sionless parameter 7 in Eqs. (3), (4). In deriving Eq. 
(15) the relations given by Eq. (7) were used. 

From this result one can easily deduce two character- 
istic features of this effect. 

1) The total resistance change between its zero-field 
value, i?(0) and its value at strong enough field (Qt s >> 
1), R(^)=R Q : 



R(oo) - R(0) _ 2 



R(0) 



= 7 



tanh A 
""A - ' 



(16). 



For a narrow sample, A << 1 the overall relative change 
of resistance is equal to 7 2 , which gives a nice way to 
determine experimentally the fundamental parameter 7. 
For wide samples the relative change is j 2 (2L s / L). 

2) The shape of the magnetoresistance curve. Wc intro- 
duce the notation p(B) for the normalized relative mag- 
netoresistance. Then 



P(B) 



R(B) - R(0) 
Woo) - R{0) 



= 1 - Re 



tanh(«A) 



k tanh A 



(17) 



For a wide sample, A >> 1, the width of the magne- 
toresistance curve is determined by the condition Qt s <~ 
1). In this case 



P(B) 



Re(-] = I - 



1 + Vl+x 2 



2(1 



1/2 



(18) 



Note that at x = Qt s » 1 the function p{B) approaches 
its maximum value very slowly, as 1 — \ j\f2x. 

Figure 2 presents numerical results for p{B) calculated 
from Eq. (17) for different values of A = L/(2L S ) to- 
gether with the curve given by Eq. (18) for A >> 1, 
which is a good approximation already for A = 1.5. 

For narrow samples (A << 1) the magnetic field de- 
pendence becomes much weaker. The reason is that 
along with the spin relaxation time r s , there is another 
characteristic time, = r s A 2 = L 2 /(4D), which is the 
time of diffusion on a distance L/2. For narrow samples 
(Td < t s ), it is this time, rather than r s , that determines 
the width of the Hanle curve, because the spin polariza- 
tion is destroyed by diffusion faster than by spin relax- 
ation. Accordingly, the width of the magnetoresistance 
curve will now correspond to fir^ ~ 1, i.e. it will be 1/A 2 
times broader compared to the case of a wide sample. 

To unify the two limiting cases, in Fig. 3 we re-plot 
p(B) as a function of the parameter fir*, where t* is 
the effective time during which the spin is destroyed be- 
cause of the combined effect of spin diffusion and spin 
relaxation: 
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FIG. 3: Normalized magnetoresistance as a function of the 
parameter £2t*, where 1/r* = l/r s + I/t^. i - in the limit 
A >> 1, 2 - in the limit A << 1. Other curves correspond 
to A = 0.4, 0.8, 1.2, 1.6, and 2. One can see that all curves 
practically coincide. 



Figure 3 shows that as a function of this parameter there 
is a quasi-universal curve, since the results for the limit- 
ing cases of a narrow and a wide samples are very close. 
Thus Eq. (18) can serve as a good interpolation formula 
for the general case, provided the variable x is replaced 
by fir*, instead of Qt s . The high- field limit is always 
approached as 1/ \/B. 

The above results for magnetoresistance are similar to 
those obtained previously [15| for the Hanle effect in the 
case when spin polarization is inhomogeneous and spin 
diffusion is important. 

From the experimental results for 3D [l(| and 2D Q 
GaAs one can estimate 7 ~ 10~ 2 . Kimura et al [l6| find 
7 = 3.7 • 1CP 3 for Platinum at room temperature, so that 
in these two cases a magnetoresistance due to spin ac- 
cumulation on the order of 10 -4 and 10 -5 , respectively, 
can be expected. The characteristic feature, allowing to 
identify this effect, is the specific form of the magnetore- 
sistance curve, as well as its strong dependence on the 
sample width, when it becomes comparable to the spin 
diffusion length. 

Because of the high precision of electrical measure- 
ments, magnetoresistance might provide a useful tool for 
studying the spin-charge interplay in semiconductors and 
metals. 
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